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Abstract 

We investigate the energetic ground states of a model two-phase system with 
1/r 3 dipolar interactions in two dimensions. The model exhibits spontaneous 
formation of two kinds of periodic domain structure. A striped domain struc- 
ture is stable near half filling, but as the area fraction is changed, a transition 
to a hexagonal lattice of almost-circular droplets occurs. The stability of the 
equilibrium striped domain structure against distortions of the boundary is 
demonstrated, and the importance of hexagonal distortions of the droplets 
is quantified. The relevance of the theory for physical surface systems with 

elastic, electrostatic, or magnetostatic 1/r 3 interactions is discussed. 
68.35.Rh, 68.35.-p, 68.10.-m, 75.30.Pd 
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I. INTRODUCTION 



There are a variety of interesting circumstances under which the spontaneous formation 
of domain structures has been observed or predicted in two-dimensional systems as a result 
of long-range electrostatic, magnetic, or elastic interactions. For example, several works 
have shown that electrostatic dipole-dipole interactions are responsible for stabilizing the 
domain patterns observed in Langmuir monolayers at the water/air interface Similarly, 
it has been suggestedHQ that work-function variations could lead to the stabilization of 
periodic domain structures on metal surfaces such as those observed for partial coverages 
of oxygen on Cu(llO).! Closely analogous is the case of thin ferromagnetic films for which 
the magnetization is normal to the plane of the film. Similar effects arise for surface 
segregation into phases which are not inherently dipolar, but which have different electric 
or magnetic susceptibilities, in an applied field. For example, dramatic effects have been 
observed in planar-confined ferrofluid/ water mixtures in magnetic fields. Finally, it has 
been understood for some time that surface stress discontinuities at domain boundaries 
can stabilize domain structures.^ This effect explains the formation of the herringbone 
reconstruction of Au(lll),il provides an alternate explanation of the Cu(110):O domain 
structuresj'i and has far-reaching consequences for Si(100) surfaces.ili^ A review of some 
of these effects appears in Ref. |25|. 

In the above cases, most of the essential physics can be captured in a model with two 
distinct phases A and B, sharp A-B boundaries having energy per unit length 7^ independent 
of orientation, and 1/r 3 interactions of strength Kd between a "dipole density" ipA or ip& in 
domains of type A or B respectively. A model of this form, or a closely related form, has been 
considered previously by many authorsBSBiiES For example, for the case of Langmuir 
films, if ipA and ipB are the screened electric dipole density normal to the water/air interface 
in the absence and presence of the Langmuir layer, respectively, then Kd is just 2e/(e + 1) 
(Gaussian units) ,@ where e is the dielectric constant of the substrate (water). For the case 
of metal surfaces, e — > 00 so that K d = 2, and Amp is just the work function.! The magnetic 
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case is analogous to the electric one. For surface stress effects, ip is just the surface stress, and 
K<i is related to a bulk elastic compliance. (In this case, ip and K^/r 3 have extra tensorial 
indices and the details become more complicated, but the scaling behavior is the same.Si!) 

Given the wide range of interesting phenomena which fall within the scope of this model, 
it is remarkable how little is known theoretically about its behavior. Even the ground-state 
structure as a function of area fraction / = /a = 1 ~ /b has not been convincingly estab- 
lished. The case of the simple striped domain structure is well understood,!l!00'l^0H but 
the stability of this domain structure to oscillatory fluctuations of the domain boundaries 
remains controversial. Marchenko has predicted that such an instability should occur, spec- 
ulating that it may lead to an ordered droplet phase.lll However, recent work of Kashuba 
and PokrovskyB and Abanov et aZ.@ does not appear consistent with such an instability. 
(The demonstration of an instability for a single isolated stripe of non-optimal widthEl does 
not necessarily imply instability of the equilibrium periodic lattice of stripes as considered 
here.) Numerical calculations of the relative energy of droplet and striped phases led to 
a predicted phase transition sequence droplet — > striped — > inverted-droplet as a function 
of area fraction /, with transitions at / ~ 0.28 and 0.72.1 These phases are illustrated in 
Fig. 1(d), (a), and (f), respectively. However, the question of the stability of the striped 
phase to sinuous distortions such as those shown in Fig. 1(b) and (c) was not addressed in 
this work. Also, the effects of deviations of the droplet shape from perfect circularity, as in 
Fig. 1(e), were not assessed. 

The purpose of this report is to establish firmly the correct equilibrium phase sequence as 
a function of area fraction within this model. We first develop the formalism to describe the 
model in Sec. II. Then in Sec. Ill, we carefully locate the transition from the striped to the 
hexagonal-droplet-lattice geometry, and study the coexistence of these two domain struc- 
tures. In Sec. IV we present numerical calculations which demonstrate that the equilibrium 
striped phase is stable against undulations of the domain walls at all area fractions, and in 
Sec. V we show that the expected hexagonal distortions of the droplets in the droplet phase 
are negligible. Thus, we confirm the basic phase sequence suggested previously! We discuss 



the limitations and implications of the work in Sec. VI, and summarize our conclusions in 
Sec. VII. 



II. DIPOLAR MODEL 



A. Formulation 



The energy in this model takes the form 



E = lb L b + ^fd 2 rf d 2 r> 




(1) 



where L b is the boundary length, and the boundary energy % is assumed to be independent 
of orientation. As noted above, is the dipolar coupling constant, and ip(r) takes constant 
values ipA or %ps when r lies in an A or B domain, respectively. (We work here at fixed area 
fraction /; one can just as easily work at fixed chemical potential /i by introducing the grand 
potential Q = E — [if A, where A is the system area and fi reflects the free energy difference 
per unit area between the two phases.) 

If desired, the second term above can be converted into a double line integral over the 
boundaries,0SEl and the energy relative to that of uniform phases A and B becomes 



Here jd = Kd{ipA — V'b) 2 , which has units of energy per unit length like 7&. When multiple 
boundaries are present, a double sum over boundaries is understood, and the sense of the 
line integrals is to be kept uniform (e.g., domain A on the right). 

For periodic domain structures, it is convenient to convert to reciprocal space, in which 
case Eq. (P or Eq. (Q) becomes^ 



We use F, as opposed to E, to represent an energy density (energy per unit area), and n b is 
the domain boundary density (length per unit area). The G are the reciprocal lattice vectors 




(2) 



AF = 76 n 6 - 7d 5>G|e(G)| 2 



(3) 



G 
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of the periodic domain structure under consideration, and 6(G) is the Fourier transform of 
the function which takes values or 1 inside domains A or B, respectively (or vice versa). 
Both 7b and 7^ are assumed to be positive, so that that the first term in Eq. (fj) or Eq. (|3]) 
suppresses domain formation, while the second promotes it. 

In each of Eqs. (0-0), logarithmic divergences occur at small length scales. These may 
be removed by introducing a microscopic cutoff a, which may have the interpretation of an 
atomic lattice spacing, a domain wall width, or a film thickness, depending on the physical 
context. It can be shown to be equivalent to introduce a into Eq. ([]]) (for sufficiently small 
a) in the form of a Lorentzian broadening (x 2 + a 2 ) -1 of the domain walls, or into Eq. (0) 
in the form of a cutoff that omits |1 — l'| < a from the integrals, or into Eq. (|3|) in the form 
of a damping factor exp(— 2Ga). 

If an appropriate value of a is not know a priori, it may be assigned an arbitrary value, 
as long as % is then fixed with respect to a known structure such as a single isolated stripe 
or disk. From Eq. (0), the energy of a single stripe of width w is, per unit length, 

dy f Y dy 



A£4 ri pe = 2 76 - 2 7d fim 



Y-KX 



a y Jo ^Jy 2 + w 2 



2 7b -2 7(i ln(^) , (4) 



,2a 

and the energy of a single circular droplet of radius R is 

™ cos 9d9 



AE disk = 27ri?7 6 - TTR-f d [ 

J a 



2ttR 



/r sin(0/2) 



The boundary energy 7^ should be regarded as being defined experimentally through Eq. 
(|) or ([D after a has been chosen. There is no real arbitrariness in practice, since a and 7^ 
enter all subsequent formulas only in the particular combination 7& + 7^1n(a). 

For later reference, it is useful to note that if the energy per unit minority-phase area 
is minimized for the isolated stripe or disk, one obtains equilibrium sizes il>o = 2l /ir and 
-Ro = e 2 /o/47r, where we have introduced a new length scale 
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l = ira exp ( — + 1 ] . (6) 

V7d / 

(As we will see in the next subsection, l has the interpretation of the stripe width in the 
equilibrium stripe phase at half filling.) The equilibrium energies per unit minority-phase 
area are 

AF stri pe = ~7T-^ (7) 

and 



AF disk = -— f . 8 
e 2 l Q 



The latter is ~8% lower in energy. Thus, circular droplets are the favored structure in the 
dilute limit. 



B. Periodic domain structures 

We have found Eq. @ to be most convenient for calculating the energy per unit area of 
periodic domain structures. In Ref. ^ it was shown that any given domain structure can be 
characterized by a dimensionless and scale-invariant constant I such that 



AE = %n b + j d n b [I + \n(nn b a)] 



(9) 



/ depends only on the geometry, and not the scale, of the domain structure, and can be 
written 



/ = lim 



- — G 1 ( G) I V 2Ga - ln(wi) 
n b G 



(10) 



The representation 



/ = lim 



7T 



]TG|9(G)|V 4G a -ln(7m 6 a)+7/2 
n b G 



(11) 
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is entirely equivalent, but has improved convergence properties (7 is Euler's constant). a Min- 
imizing Eq. (H) with respect to n b , and making use of Eq. (^), one finds that the equilibrium 
value of rib is given by 

n^ = l exp(I) , (12) 

at which the energy density is 

AF = -^exp(-J) . (13) 

Thus, regardless of domain geometry, the uniform phase is always unstable to the formation 
of a domain structure, although the length scale n^ 1 of this domain structure depends 
exponentially on the parameters of the model through Eqs. @ and (|13|). Note that ■jd/lo m 
Eq. (|13|) sets the energy scale, while loAF/jd = — exp(— J) plays the role of a dimensionless 
energy function. 

It is clear from Eqs. ( |12"D and (|3~3|) that the domain structure which minimizes I will have 
the lowest energy, and will thus be the physical equilibrium structure at zero temperature. 
For a simple striped phase characterized by area fraction /, as illustrated in Fig. 1(a), one 
finds 1(f) = — In sin(7r/).!il This has an absolute minimum / = at / = 1/2 for stripes of 
equal width; for this case n^ 1 is just lo, showing that Iq is just the equilibrium stripe width. 

The principle problem to be addressed next is whether other periodic domain structures, 
such as wavy stripes or droplets, might have lower / than the simple striped structure, at 
least in some range of /. 

III. COMPETITION BETWEEN STRIPE AND DROPLET PATTERNS 

In Ref. |6], it was shown that a hexagonal lattice of circles of one phase on a background of 
the other phase, shown in Fig. 1(d), becomes energetically favorable relative to the striped 
structure in the vicinity of / < 0.28. (Of course, by symmetry, the inverted droplet phase 
of Fig. 1(f) is then favored for / > 0.72.) We have repeated the calculation of 1(f) for the 
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droplet lattice, this time using the analytic representation of the 2D Fourier transform of a 
disk 

9(g) = ^MqR) (14) 

( Ji is a Bessel function). The results are shown in Fig. |2|, where the scale-optimized energies 
of Eq. fll3l ) of striped and droplet structures are plotted in dimensionless form vs. area 
fraction /. It can be seen that the energies of the droplet and striped domain patterns are 
very close over most of the range < / < 1/2, the deviation at / = 1/2 being no more that 
~5%. The crossover between the two curves is found to occur at f c = 0.286. The slopes of 
the curves at / — > correspond to the critical values of chemical potential \x at which the 
minority phase disappears. These are given by —n = —3.142 and — 87r/e 2 = —3.401 as per 
Eqs. (0) and (||) for striped and droplet domains, respectively, and differ by only ~8%. 

Very near the critical area fraction f c , the system can lower its energy by phase sepa- 
rating into superdomains of striped and droplet phases having / slightly greater and less 
than f c , respectively. The coexistence region is determined by the usual common tangent 
construction applied to the energy function plotted in Fig. |2|. The result of this construction 
is shown in the inset, where an irrelevant linear function has been added to aid visibility. We 
find that the coexistence region is delineated by a droplet phase at f c \ = 0.273 in equilibrium 
with a striped phase at f c2 = 0.299 at fi c = —1.855. 

In the region of phase separation, f c i<f< f C 2, the superdomains themselves should in 
principle order into a periodic domain superstructure, by virtue of the fact that they have 
slightly different dipole densities ipes = V'A + /(V'b — V'a)- The theory of Sec. |I A| applies, 
now with 7^ = Kd(ipB — 4>A) 2 Af 2 , where Af = / c2 — / c i=0.026, and is interpreted as the 
energy per unit length of a superdomain boundary. Thus a striped superstructure is to be 
expected near the middle of the phase separation region, and droplet superstructures might 
occur near f cl or f C 2^ (If so, it is even conceivable that these super-phases could phase 
separate near coexistence into periodic super-super-structures, etc.) However, the energy 
scale for the superordering would be extremely weak (Af 2 < 10~ 3 ), and it seems doubtful 
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whether these super domain effects could be observed experimentally. 



IV. STABILITY OF STRIPED PHASES 



We now wish to consider the stability of the striped domain phase to perturbations of the 
boundaries. We will consider both in-phase and out-of-phase fluctuations of the boundary, 
as sketched in Figs. 1(c) and 1(d), respectively. We let the stripes be propagating along the 
x-direction, with unit periodic repeat distance in the ^-direction, and a boundary deviation 
of the form of 



The problem is to compute I(f, k, A), where / is the area fraction, and A is the amplitude 
and k the wavevector of variation. By symmetry, I(f, k, A) can be expanded in even powers 



This equation defines the stiffness a(f, k) of the domain wall with respect to sinuous dis- 
placements. If a negative value of a were found, it would indicate instability of the striped 
domain structure. 

We first develop a method for calculating a(f,k) directly. The perturbed structure is 
periodic with lattice vectors 2nxjk and y, and reciprocal lattice vectors G = kmx + 2imy, 
so that the Fourier components appearing in Eq. (^) can be written 0(G) = @(m, n). 
Expanding in powers of A and keeping only terms up to order A 2 , we find for the case of 
in-phase variation of the domain boundaries 



y = A sin(kx) 



(15) 



of A: 




(16) 




(17) 



nn 



and 




(18) 
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with all other Fourier components vanishing. Equation ( |T0"D then leads to 



a 



(/,fc) = E sin W) 



n=l 



-4mra i 



2n 



+ A-K 2 n - 2irak 2 ) 



2nJ(2Trn) 2 + k 2 e ~W( 2 ™) 2 + k2 a 



(19) 



Introducing k = k/2ir and taking the limit a —>■ 0, 



«(/> K ) = E 471-2 sin 2 (7rn/) 



n=l 



h n — y/n 2 + k 2 

2n 



(20) 



which can be proved by careful expansion of the exponential factors in Eq. fll9f) . The three 
terms appearing in the bracket above result from the increase of rib, the decrease of the 
0(0, n) components, and the increase of the 0(±1, n) components, respectively, with A. 
Since the first two terms are always larger in magnitude than the third term, we can see 
that a is always positive. The results for ol{k) are plotted for / = 1/2 in Fig. |3|. In the limit 
k — > 0, ot is directly proportional to k 4 , while for very large k, a scales as k 2 \ii(k). These 
limiting behaviors are illustrated in Fig. 0. In particular, we find that a ~ 1.0527t 2 k 4 for 



0, and a ~ 7r k \u(k) for k large, at / = 1/2. 



We use a'(f, k) to denote the corresponding stiffness for the case of out-of-phase variation, 
Fig. 1(c). For this case, the m = elements of are still given by Eq. flTT|), while 
becomes 



0(±l,n) = ie- mn Acos(7rnf) 



(21) 



We now find, as a — > 0, 



a 



n=l 



— COS' 



K 



sin 2 (7rn/) ( \-n 

\2n 



(7mf)Vn 2 + k 2 



(22) 



The series is conditionally convergent and has to be summed carefully. The results for 
/ = 1/2 are shown as the dashed line in Fig. |3]. It can be seen that a' > a (and thus a' > 0) 
over the entire range of wavevector; this was also found to be the case at other area fractions 
/. For / = 1/2 the minimum occurs at about k = 47r/3. 
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Very recent work of Abanov et a/.lll contains a similar stability analysis which appears to 
confirm these results for the special case of half filling. Our Eqs. (|20| ) and ( |22|) with / = 1/2 
correspond to their Eq. (A23) with p x = and tt/L, respectively. 

We have further checked the above results numerically by carrying out a direct numerical 
evaluation of Eq. (|TT|). For this we do not want to rely on an analytical Fourier transform 
of O(r), so we begin by specifying 0(r) directly on a real-space mesh covering the unit cell 
of interest. To avoid artifacts of the discreteness of the lattice, the function G'(r) that we 
actually evaluate on the mesh is a Gaussian-smeared version of O(r), 



§(r) = i 



erf 



(23) 



,2V2a, 

where erf (x) is the error function and \d\ is the distance from the domain boundary to r, the 
sign of d depending on whether r is in domain A or B. Since this corresponds to convoluting 
with a Gaussian function in real space, we have in reciprocal space 

9(G) = e(G)e" 2G2a2 . (24) 

In practice, 0(G) is computed by fast Fourier transform. Then the sum in Eq. ([□]) can 
be computed simply as J2g ^l@(G)| 2 . We found that good numerical accuracy was ob- 
tained using ~1024-2048 grid points in the ^-direction, and an equal mesh density in the 
x-direction. The calculated values of 1(A) were then fitted to a quartic polynomial in A 
in order to determine a or a'. The values determined by this approach match very well 
with those computed from Eqs. (|20|) and (0), as shown in (Fig. |3|) for / = 1/2. The errors 
are typically ~3% or less. Moreover, by using this approach, we confirm that the surface 
energy change per unit area increases monotonically with A in all cases considered. This 
provides additional evidence, above and beyond the linear-response 0(A 2 ) analysis, that 
the equilibrium in-phase and out-of-phase striped domain structures are always stable with 
respect to domain boundary variations. 
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V. HEXAGONAL DISTORTION OF DROPLETS 



We now consider in more detail the domain phase consisting of a hexagonal lattice of 
droplets. There is no reason to expect the droplets in this phase to be exactly circular; 
the interactions with neighboring droplets should give rise to some amount of hexagonal 
distortion, not considered in the previous section, as shown in Fig. 1(e). Are these distortions 
important? For example, do they lower the energy of the droplet phase significantly with 
respect to the striped phase, perhaps even eliminating the striped phase from the phase 
diagram altogether? On the contrary, we show here that the effects of the distortion are 
insignificant. 

We thus consider boundary distortions of the form 

r = r + A cos(6#) , (25) 

where r and 9 are polar coordinated measured from the center of the droplet, r is its 
unperturbed radius, and A is the amplitude of variation. Unlike the case of the striped 
phase, there is no symmetry which requires /(/, A) to be an even function of A, so we 
expect its Taylor expansion in A will contain linear as well as higher-order terms. Because 
there is no simple analytic method to calculate ©(G) for given / and A, we again specify 
@(r) first on a real space mesh, fast Fourier transform to reciprocal space, and evaluate 
Eq. ([0]) in G-space. We find that a 2048x2048 mesh in the parallelogram unit cell gives 
sufficient numerical accuracy to calculate /(/, A). Linear and quadratic coefficients in A are 
obtained from a simple fit to a series of calculations using different values of A, and used to 
locate the equilibrium distortion A min and the corresponding energy change AI = J min — I 
at v4 m i n . 

We have carried out such calculations at many values of /, but it is sufficient to report the 
results at / = 1/2 (where the striped phase has maximum stability) and at / = 0.28 (near 
the crossover between the droplet and striped phase). For the cases of / = 0.5 and / = 0.28, 
we find Anm/r = 3 x 10~ 3 and 2 x 10" 3 , and AI = -1.3 x 10~ 4 and -2.6 x HT 4 , 
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respectively. We find A m - m < for all /, which with our conventions indicates that the 
boundaries try to avoid each other. This causes the portions of the droplet boundary which 
directly face other droplets to become flattened, as indicated (in exaggerated fashion) in 
Fig. 1(e). In principle, higher harmonics such as cos(12#) should be considered in Eq. (p5|). 
but given the smallness of the fundamental 66 distortion, these are virtually certain to be 
negligible as well. 

The effect of this AI on the relative stability of the striped and droplet domain phases 
is miniscule. We estimate that the critical filling fraction f c which separates the droplet and 
striped phases is shifted by only ~ 10~ 4 due to the distortion of the droplets. Thus, the 



exact values of f c , f c \ and f C 2 of Sec. |T| are only very slightly shifted, and the description 
of the transition is virtually unchanged. 



VI. DISCUSSION 

While the model that we have studied is motivated by the experimental work summarized 
in the Introduction, it should be emphasized that extensions of various kinds would probably 
be needed to make real contact with experiment in most cases. 

Perhaps the closest experimental realization of the model can be found in the case of 
Langmuir layers at the air-water interface. Experimental work on this system has shown 
indications of the presence of both droplet and striped phases, although it is not clear that 
the transition between them is connected experimentally with a change in area fractionUHl 

Probably the most serious limitation of the model is the assumption that the boundary 
energy 7^ is independent of orientation. While this is clearly correct for the case of Langmuir 
layers at the air-water interface, and for ferrofluid mixtures. But it is clearly not correct 
at crystal surfaces. For example, the striped domains observed on the Cu(110):O surfaces 
always run along the [001] direction,i presumably because [110] boundaries are much lower 
in energy than [110] ones. Anisotropy in 7^ will always favor striped domains over droplet 
or inverted-droplet domains, so the effect of a weak anisotropy will be to shift the critical 
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area fractions for droplets closer to / = and / = 1. A strong anisotropy will eliminate the 
droplet and inverted-droplet phases from the phase diagram altogether. 

Another limitation is that in some of the physical situations discussed, the idealization 
of a pure two-dimensional model is not quite appropriate. For example, for the work on 
confined ferrofluids and epitaxial ferromagnetic films, the film thickness plays a definite 
role, giving rise to specific forms of the effective interaction between A and B domains and 
eliminating the need to introduce an ad hoc real-space cutoff a. In the case of Coulomb 
interactions in an electrolytic background,^ screening can convert the r -3 interaction into a 
short-ranged form, while for the elastic case, the surface stress ip and the elastic interaction 
Kd/r 3 carry extra tensorial indices. 

Also, we have assumed that the dipole density ip always takes on two discrete values 
ipA and ipB- In some cases we may prefer to think in terms of a model in which there is 
an intrinsic energy density F(ip) having local minima at ip — ipA and ips, but for which 
variations of ip from these minima are allowed. For the case of Langmuir layers at the 
air-water interface, for example, the molecular density in the Langmuir layer is sure to be 
somewhat compressible. 

Moreover, in some cases our characterization of the domain boundary as infinitely sharp 
may be inappropriate. For example, in the case where phases A and B arise from phase 
segregation, the boundary is expected to become broad near the phase-separation critical 
temperature. In such cases also, it is appropriate to treat ^(r) as a continuous field, using a 
Landau-type expansion to describe the energy of the system and the appearance of domain 
structures.!^ 

For atomic-scale magnetic systems, it might be appropriate to go over to a descrip- 
tion in terms of lattice spin models, instead of insisting on a continuous r-space. The 
simplest such model would be an Ising model on a square 2D lattice with ferromagnetic 
nearest-neighbor and antiferromagnetic r~ 3 long-range exchange interactions. (Unless a 
specially-tuned second-neighbor interaction is introduced into such a model, the anisotropy 
of the boundary energy is again likely to favor striped phases aligned along the Cartesian 
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directions.) We are not aware of previous studies of such a model. 

Even assuming the formulation of the model is correct for the system of interest, it would 
be very desirable to generalize the above zero-temperature theory to finite temperature. The 
energy scale for defects in the domain structures (e.g., vacancies or interstitials in the droplet 
phase, stripe terminations in the striped phase) is 7^0, so it is natural to introduce a dimen- 
sionless temperature t = ksT/jJo- Thus, the problem reduces to mapping out the phase 
diagram in the / — t (area-fraction-temperature) plane. This remains an ambitious program 
for future work. However, we can make the following speculations. We certainly expect 
a solid-to-liquid melting transition in the droplet phases, quite possibly via the Kosterlitz- 
Thouless mechanism, with a narrow hexatic phase interposed. The melting temperatures 
should vanish at / — > or 1, since the interactions between droplets becomes very weak 
in those limits. Theoretical considerations suggest that the striped domain structure will 
melt at any non-zero temperature, giving rise to a 2D nematic phase characterized by an 
orientational correlation length and having exponential and algebraic decay of positional 
and orientational order, respectivelyiEii (Further discussion for the case of anisotropic 75 
appears in Refs. |T5] , ^B[ ) The nature of the low-temperature part of the phase diagram near 
f c is especially unclear. 

Finally, even assuming that the model is correct and that we know its equilibrium behav- 
ior, there may be many circumstances under which the kinetic behavior is more important 
that the thermodynamics. In fact, we may expect sizable energy barriers to the coalescence 
or pinching off of droplets, or the formation and annihilation of stripe crosslinks. Some 
studies of the dynamics of phase separationEE and fingeringEiiB in dipolar systems have 
already appeared. Nevertheless, even in cases where kinetic effects dominate, it may be 
of interest to understand the equilibrium phase diagram first, in order to understand the 
driving forces for the dynamic effects observed. 
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VII. SUMMARY AND CONCLUSIONS 



In summary, we have carried out a theoretical study of the stability of periodic striped and 
droplet domain structures for a simple two-dimensional dipolar model at zero temperature. 
It is confirmed that the striped structure is stable near area fraction / = 0.5, with transitions 
to droplet and inverted-droplet structures at / = 0.286 and / = 0.714, respectively. A small 
phase-separation region near the striped-to-droplet transition was identified. The stability of 
the striped structure against sinuous displacements of the domain boundaries was confirmed, 
and the hexagonal distortion of the droplets in the droplet phase was quantified and found 
to be negligible. 

The characterization of the finite-temperature phase diagram of the model, and its kinetic 
evolution, remain as challenging problems for future study. Furthermore, application to 
some experimental systems may require refinements of the model, such as an orientation- 
dependence of the domain-wall energy, a finite film thickness, a finite domain-wall thickness, 
a compressibility of the dipole density, or addition of tensorial indices for the elastic case. 
Nevertheless, it is hoped that the present work will provide a firm foundation for future 
developments. 
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FIGURES 

FIG. 1. Striped domain phase: (a) unperturbed, (b) with in-phase, and (c) with out-of-phase 
boundary displacement, Eq. (|T^). Droplet phase: (d) unperturbed; (e) with hexagonal boundary 
displacement, Eq. (|25|). (f) Inverted droplet phase. 

FIG. 2. Comparison of dimensionless energy per unit area, optimized with respect to scale, 
for striped (solid line), droplet (long dashed), and inverted droplet (short dashed line) domain 
structures, as a function of area fraction /. Inset illustrates common tangent construction which 
determines the phase coexistence region. 

FIG. 3. Stiffness a vs. reduced wavevector k for striped phase at filling / = 1/2. Solid and 
dashed lines denote in-phase and out-of-phase boundary variations, as computed from Eqs. (|20| ) 
and ([22]), respectively. Triangles and circles indicate the corresponding results obtained numerically 
by finite differences from calculated values of /(/, k, A) using the method of Eqs. (^3|)-([2"4|). 

FIG. 4. Stiffness vs. wavevector for in-phase boundary variations, from Eq. (|20|), illustrating 
asymptotic form for large and small k. Results are expressed as a/n 2 vs. k 2 for three values of 
area fraction /. 
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